Wavelength dependence of high-order harmonic generation with the ionization and 
ponderomotive energy controlled by an extreme ultraviolet pulse 
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We theoretically study the scaling with the driving wavelength A of the high-order harmonic gener- 
ation (HHG) under the simultaneous irradiation of an extreme ultraviolet (XUV) pulse. Surprisingly, 
when the cutoff energy and ionization yield are fixed, the harmonic yield is nearly independent of 
A. We identify its origin as the combination of the initial spatial width of the states excited by the 
XUV pulse, making the wavepacket spreading less prominent, and the shallowing of the ionization 
potential, which suggests complex nature of the wavelength dependence of HHG. 

PACS numbers: 42.65.Ky, 32.80.Rm, 42.50.Hz, 32.80.Fb 
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High-order harmonic generation (HHG) represents one 
of the best methods to produce ultrashort coherent light 
covering a wavelength range from the vacuum ultra- 
violet to the soft X-ray region. The maximal har- 
monic photon energy Ec is given by the cutoff law 
Ec — Ip -\- 3.17C/p m, where Ip is the ionization po- 
tential of the target atom, and Up[eW] ~ F'^/Alo'^ — 
9.337 X IQ-'^^ I [W/cm^] (A [mhi])^ the ponderomotive 
energy, with F , I and A being the strength, intensity and 
wavelength of the driving field, respectively. Since Up 
scales as A^, a promising route to generate harmonics of 
higher photon energy is to use a driving laser of a longer 
wavelength. Thus, the laser wavelength (A) is an effective 
control knob for the ponderomotive energy, and hence the 
cutoff. This has motivated HHG experiments with high- 
power mid-infrared (MIR) lasers (20,1^]. Using a 1.55 
/Ltm driving laser field from an optical parametric ampli- 
fier for example, Takahashi et al. P have recently 
succeeded in generating harmonics with a photon energy 
of 300 eV from Ne and 450 eV from He gas, which lie 
well in the water-window region. 

Under such a circumstance, the dependence of the 
HHG yield on A has become an issue of increasing in- 
terest. Although it had been commonly assumed that 
the HHG efficiency scaled as A~'^ due to the spreading of 
the returning wavepacket 0, recent theoretical 0, U H, 
lol . 11 1 as well as experimental^,!!] studies have revealed 
much stronger dependence of esc A~^ with 5 < a: < 6, 
which would significantly reduce the HHG yield by MIR 
lasers. It is considered that the additional factor A~^ is 
of an apparent nature stemming from the distribution of 
the HHG energy up to the cutoff (cx A^) [!,[§], though 
the precise physical origin of the scaling law has not been 
fully understood yet. 

While most of the experiments are conducted with a 
driving laser of a single wavelength, the control of HHG 

using XUV pulses has also been discussed H E, m. 



For example, Schafer et al. [I4I showed that the delay 
of attosecond pulse trains can be used to microscopically 
select a single quantum path contribution. On the other 
hand, Ishikawa [l3, 14| has theoretically shown that the 
irradiation of the XUV pulse with a photon energy fiujx 
smaller than Ip can boost the ionization Yj and harmonic 
yield by orders of magnitude; the XUV pulse facil- 
itates optical-field ionization by promoting a transition 
to (real or virtual) excited states. This effect has been 
experimentally demonstrated by the use of mixed gases 
[15| . and its application to single attosecond pulse gen- 
eration has been proposed [3]. One of the remarkable 
features of this effect is that Yj increases in proportion to 
the XUV intensity without affecting the cutoff energy de- 
termined by the driving infrared pulse (Fig.[T]). Thus, the 
addition of XUV pulses can be viewed as a tool to enable 
independent control of A, E^, and Yj: for a given value 
of A, Ec can be adjusted through / and then Yj through 
the XUV intensity. This would provide the investigation 
of the A-dependence of the HHG with additional degrees 
of freedom. 

The above consideration has motivated us to theoret- 
ically investigate the driving-wavelength-dependence of 
HHG with the XUV control of Ec{Up) and Yj. For the 
case of the driving laser pulse alone, if we fix Ec at each 
driving wavelength, the driving intensity is lowered with 
an increasing wavelength, leading to the drop of Yj, which 
in turn largely affects the HHG efficiency. The addition 
of an XUV pulse of appropriate intensity, however, can 
adjust y to a constant value, and, then, we would expect 
oc A"'^ scaling due to the wavepacket spreading. Our re- 
sults based on numerical solution of the time-dependent 
Schrodinger equation (TDSE), however, show that the 
harmonic yield is nearly independent of A at fixed pon- 
deromotive energy and ionization. Using the Lewenstein 
model we identify the origin of this surprising fea- 
ture as the combination of the initial spatial width of the 



2 



wave function and shallowing of the effective ionization 
potential, indicating complex nature of the A-dependence 
of HHG. 

To study the single-atom response under a combined 
driving laser and XUV pulse, we solve the TDSE in the 
length gauge. 
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Vir) + z[E{t)+Ex{t)] 



(1) 

for a model atom in the single active electron approxi- 
mation, represented by an effective potential 17|, 



V{r) = -[l-f ae-'' + {Z - 1 - a)e-'^'']/r, (2) 
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— Fund. (LSxIO^^W/cm^) & XUV (¥,=0.31%) 
Fund, only (5x1 o'^W/cm^) (Y|=0.31%) 

< (1.6x10"W/cm^) (Y|=1.7x10"Vo) 




100 



2 3 

(^-ipyup 

200 300 
Photon enegy (eV) 



4 



400 



500eV 



where Z denotes the atomic number. For He, we use 
parameters Z — 2, a — 0, and /3 — 2.157, which faith- 
fully reproduce the eigenenergies of the ground and the 
first excited states. E{t) = Ff{t)smLut is the driving 
optical field, with F being the peak amplitude and f{t) 
the envelope function corresponding the Gaussian pro- 
file with a full width at half maximum (FWHM) of 35 
fs. Exit) = Ffx{t)sinujxt is the XUV field, with Fx 
being the peak amplitude. The harmonic spectrum is 
calculated by Fourier transforming the dipole accelera- 
tion, and the HHG yield is defined as energy radiated 
from the target atom per unit time [l8| integrated for a 
fixed range of photon energy hujh, specifically from 30 to 
60 eV. 

Figure [T] shows the harmonic spectra from He for 
A = 1600 nm with and without the XUV field {huix — 
17.05 eV). For the case of the driving laser alone with a 
peak intensity / of 1.6x 10^** W/cni^ (blue dashed curve), 
the ionization yield Yj is very low (1.7 x 10~^%). We 
can increase Yj in two ways. First, if we augment / to 
5 X 10^'*W/cm^ (black dotted curve), Yj reaches 0.31% 
and, accordingly, the harmonic yield becomes higher, 
which is accompanied by the increase of the cutoff energy. 
Alternatively, the same ionization yield can be achieved 
by the addition of the XUV pulse with an appropriate 
intensity of 2.3 x lO^^W/cm^. In this case (red solid 
line), the cutoff remains nearly unchanged. Hence, as al- 
ready mentioned, the combination of the laser and XUV 
pulses can be used as a tool to adjust A, Ec{Up), and Y 
independently. It should also be noted that the resulting 
harmonics have an even higher yield than those from a 
driving laser of higher intensity alone (black dotted line) 
between 30 and 60 eV. The ratio between the two cases 
in this energy range is « 3.2, which is comparable with 
the ratio of Up. In addition, as is shown in the inset, 
the harmonic yield is distributed in a similar manner be- 
tween hujh = Ip and Ec, in spite of the large difference in 
driving intensity and cutoff energy. These observations 
are consistent with the idea that the additional scaling 
cx is an apparent effect due to the harmonic energy 
distribution up to the cutoff. 



FIG. 1: (color online) Upper solid curve: harmonic spectrum 
from He exposed to a 35 fs Gaussian combined driving and 
XUV pulse {htox = 17.05 eV), the former (A = 1600 nm) 
with a peak intensity of 1.6 x 10^''W/cm^ and the latter 
2.3 X lO^^W/cm^. Middle dotted and lower dashed curves: 
harmonic spectra for the cases of the driving pulse alone, with 
an intensity of 5 x 10^"' W/cm^ and 1.6 x lO'^'* W/cm^, respec- 
tively. Inset: replots of the upper two curves in terms of 

{hujh - ip)/Up. 
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FIG. 2: (Color onHne) Wavelength dependence of the TDSE- 
calculated DE harmonic yield from He between 30 and 
60 eV, for different values of Tuajx- I ~ 1.6 x 10^"* x 
((800 nm)/A)^ W/cm^ so that Up may remain unchanged, 
except for the triangles and the fitting line, for which I is 
fixed at 1.6 x 10^* W/cm^. The XUV intensity is adjusted in 
such a way that the ionization yield is 1%, irrespective of A. 

Encouraged by these results, let us now explore how 
the harmonic yield variates with the driving wavelength 
when Up{(x iX^) and ionization are kept constant si- 
multaneously by the addition of an XUV pulse. Many 
features of HHG can be intuitively and even quantita- 
tively explained by the semi-classical three-step model 
[E [3, 113 ■ According to this model, an electron is lifted 
to the continuum at the nuclear position with no kinetic 
energy (ionization), the subsequent motion is governed 
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classically by an oscillating electric field {propagation), 
and a harmonic is emitted upon recombination. The last 
step is independent of A as far as a given harmonic pho- 
ton energy range is concerned. The first step is fixed. 
Concerning the propagation step, if we neglect Ip in the 
saddle-point equations |2l[ , or equivalently, if we con- 
sider a classical motion of electron in an oscillating elec- 
tric field starting from the origin with a vanishing initial 
velocity, the phases of the field upon ionization 0j = Luti 
and recombination (jjr = Lut^ {ti,tr: time of ionization 
and recombination, respectively), characterizing quan- 
tum trajectories, are a function of hcuh/Up, hence com- 
mon for any value of A, since Up is fixed. Thus, we might 
expect that the comparison under the condition of fixed 
ionization and Up extracts the effect of the wavepacket 
spreading. 

In Fig. [2] we show the dependence of the XUV-assisted 
harmonic yield on the driving wavelength from 800 nm 
to 1.6 ^m for several different values of XUV photon 
energy hujx, including 20.964 eV resonant with the tran- 
sition to the first excited state. The peak intensity / is 
1.6 X 10^'* W/cm^ at A = 800 nm and varied so that Up 
(oc IX^) remains unchanged. The XUV intensity is ad- 
justed to yield Y = 1%, irrespective of A. We can see 
that, apart from fluctuations due to quantum-path in- 
terference [1, [^, [13, 11 , 2^ , the harmonic yield is nearly 
independent of driving wavelength, in great contrast to 
the common anticipation that the wavepacket spreading 
has a contribution cx A^^. In this figure is also shown the 
result for the driving intensity fixed at 1.6 x lO^'' W/cm^; 
the ionization yield is adjusted again to 1% through the 
XUV intensity, though it scarcely depends on A. In this 



case, reflecting the apparent harmonic energy distribu- 
tion effect, the HHG yield scales as A~^, which is much 
gentler than the usual A~^ dependence for the case of the 
driving pulse alone. 

In order to clarify the origin of this surprising fea- 
ture, let us re-examine the wavepacket spreading dur- 
ing the propagation process. The enhancement mecha- 
nisms under simultaneous irradiation of the XUV pulse 
are harmonic generation from a coherent superposition 
of states and two-color frequency mixing_^ (tunneling ion- 
ization from a virtual excited state) [ij, [l3]- The ex- 
cited states are spatially much more extended than the 
ground state. Our discussion so far as well as the com- 
mon discussion on the wavelength dependence, however, 
neglects the initial spatial width of the wave function. 
The latter can be explicitly accounted for in the Lewen- 
stein model Q if we approximate the ground state by 
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a Gaussian wave function ^(r) = (ttA^) 
where A Ip^) is the spatial width. An appealing 
point of this Gaussian model is that one can analytically 
evaluate the integral with respect to momentum in the 
formula for the dipole moment (Eq. (8) of Ref. 0). The 
spreading factor (2A^ + ir)^^/"^ in the resulting formula 
(Eq. (22) of Ref. 6]) includes the effect of the width of 
the initial state. 

Let us here extend the above discussion to the HHG 
from the superposition of the ground and an excited 
states, relevant to the enhancement mechanism 13, l3|- 
Then, following Ref. [2^, we obtain the formula for the 
dipole moment d(t) as. 



d[t) = i(AgAe)-'/2 / {2C{t,t')f/^E{t'){A{t)A{t')+C(t,t')[l-D{t,t')iA{t)+A(t'))]+C\t,t')D^{t,t')} 



exp -~-i[{Ipt - let') 



A^{t)^l + A^{t')^l - C(t, t')D^(t, t') \ 
B{t,t')] ^—^ , (3) 



r 



where Jg denotes the ionization potential of the excited 
level, Ag and Ag the spatial width of the ground and 
excited states, respectively, A{t) the vector potential, and 

B{t,t')^^j'dt"A^t"), (4) 
Cit,t') = {Al + Al+^it-t'))-\ (5) 



D{t, t') = A{t)Al + A{t')Al + i J dt"A{t"). (6) 

The factor C^/'^(t,t') describes the leading contribution 
from the wavepacket spreading. For the first excited state 



of He (/e = 3.6 eV), for example, A^ is several tens 
of a.u., hence comparable with the excursion time r — 
t — t' . This, making the wavepacket spreading relatively 
less prominent, is expected to influence the wavelength 
scaling. 

It should be noted that if we resorted to the saddle- 
point analysis (SPA) [3, [2l[ instead of the momentum 
integration, the ionization time t' would contain an imag- 
inary part Imt' w \/2Ie/ E{t') stemming from the tunnel- 
ing process. Then the spreading factor would rather read 



as A^ 



, — g , y2Ie/E{t') + ir) containing an ad- 
ditional term that can be interpreted as the width at the 
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TABLE I; Exponent x of the wavelength scaUng oc for 
various combinations of the initial spatial width Ae and the 
effective ionization potential Je. 
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tunnel exit 2^ . This tunneling contribution is automati- 
cally accounted for in Eq. Since < ^/2]^/E{t') < 
y/2Ip/ E{t') < Ag in general, the width of the excited 
state has the largest contribution in the XUV-assisted 
HHG while the initial width is neghgible for the case of 
the ground-state atom. 

The form of Eq. ([3]) suggests that the dependence of 
the HHG yield on the initial spatial width of the wave 
function and the ionization potential is rather complex. 
While these two are correlated to each other in the real 
atom, here we treat them as independent parameters and 
list in Table |l] the exponent of the power-law scaling for 
different combinations of Ag and /g, calculated with Eq. 
([3]). /e as a free parameter may be interpreted as the 
effective ionization potential defined by Ig = Ip — hux- 
It should be noted that the peak intensity is fixed at 
1.6 X 10^^ W/cm^, so the results are to be compared with 
the triangles in Fig. [2l Both larger initial spatial width 
and shallower effective ionization potential decrease the 
exponent, and their synergy leads to the surprising gentle 
wavelength scaling. 

In summary, we have investigated the driving- 
wavelength dependence of HHG under the simultaneous 
irradiation of a non-ionizing XUV pulse. The XUV pulse 
serves as a tool to provide additional degrees of freedom 
to the study of the A-dependence of HHG, with its abil- 
ity to adjust the cutoff and ionization yield independently 
and control the initial spatial width of the wave function 
and the effective ionization potential. We have shown 
that the XUV-assisted harmonic yield scales with A much 
more weakly than for the case of the driving laser alone; 
fixed Up and Y, especially, lead to a very small A depen- 
dence. According to our analysis based on the Gaussian 
model, the combination of the large spatial width of the 
states excited by the XUV pulse making the effect of the 
wavepacket spreading less prominent and the shallowing 
of the effective ionization potential is responsible for this 
unexpected feature. While both effects are described in 
Eq. ([3]) in principle, clear-cut explanation why the latter 
contributes to the gentle scaling is not at hand. The re- 
sults of the present study indicate that the A-scaling of 
HHG is not simply governed by the wavepacket spread- 
ing (cx A^'^) and the apparent energy distribution effect 



(cx A~^), but exhibits richer and more complex behav- 
ior than previously considered. There are indeed further 
open questions such as why higher-order returning tra- 
jectories have so important contribution BHBIio'l and 
whether the wavepacket spreading should give a factor 
oc A^^ rather than cx A^'^ since the spreading in the di- 
rection of the quiver motion is swept over the parant ion 
upon recoUision. Further study will be necessary to an- 
swer these questions. 
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